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We present the first proteomic analysis on the cellular 
response to severe acute respiratory syndrome-associ- 
ated coronavirus (SARS-CoV) infection. The differential 
proteomes of Vero E6 cells with and without infection of 
the SARS-CoV were resolved and quantitated with two- 
dimensional differential gel electrophoresis followed by 
ESI-MS/MS identification. Moreover isotope-coded affin- 
ity tag technology coupled with two-dimensional LC- 
MS/MS were also applied to the differential proteins of 
infected cells. By combining these two complementary 
strategies, 355 unique proteins were identified and quan- 
titated with 186 of them differentially expressed (at least 
1.5-fold quantitative alteration) between infected and un- 
infected Vero E6 cells. The implication for cellular re- 
sponses to virus infection was analyzed in depth accord- 
ing to the proteomic results. Thus, the present work 
provides large scale protein-related information to inves- 
tigate the mechanism of SARS-CoV infection and patho- 
genesis. Molecular & Cellular Proteomics 4:902-913, 
2005. 


A new type of coronavirus was reported as the causal agent 
of severe acute respiratory syndrome (SARS)' in April 2003, 
and the genome of the SARS-CoV was sequenced by several 
groups (1-3). The properties of the SARS-CoV genome was 
analyzed in depth by bioinformatic tools (4, 5). In addition, 
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several important works on the proteins of SARS-CoV have 
been reported recently, including the identification of SARS- 
CoV 3C-like protease structure (6) and the identification of 
angiotensin-converting enzyme 2 as a functional receptor for 
the spike protein (7). In our recent work, we identified all of the 
predicted SARS-CoV structural proteins, nucleocapsid (N), 
membrane (M), spike (S), and envelope (E), using proteomic 
approaches and found a novel protein, SARS-CoV 3a (8, 9). 

To uncover the mechanisms of cellular responses to the 
virus infection and identify potential drug targets of antiviral 
treatment, it is very useful to study the molecular profiling of 
virus-infected cells with high throughput and quantitative ap- 
proaches. Analysis of gene expression profiles during viral 
infection is one of the powerful approaches to probe potential 
cellular genes involved in viral infection and pathogenesis (10). 
The recent development of proteomic analytic technology 
such as differential gel electrophoresis (DIGE) (11, 12) and 
ICAT (13, 14) also provides new tools for such studies. 

As a method based on two-dimensional (2D) electrophore- 
sis, DIGE allows two or three independent samples labeled 
with different fluorescent dyes such as cyanine-2 (Cy2), cya- 
nine-3 (Cy3), and cyanine-5 (Cy5) to be run in one gel simul- 
taneously and viewed individually using the different fluores- 
cent properties of Cy2, Cy3, and Cy5, circumventing some of 
the reproducibility problems associated with 2D electrophore- 
sis and providing more accurate quantitative information 
compared with other staining methods such as silver staining 
with the dynamic range over 3—4 orders of magnitude (11, 12). 
More recently, the combination of stable ICAT, LC, and 
MS/MS has emerged as an alternative gel-free quantitative 
proteomic technology (13, 14). In ICAT analysis, two pools of 
proteins are labeled respectively with isotopically light and 
heavy ICAT regents, which are chemically identical and there- 
fore serve as a good internal standard for accurate quantifi- 
cation. Although LC-MS provides quantitative information 
based upon the relative abundances of the heavy and light 
peptides, LC-MS/MS provides qualitative information based 
upon the peptide molecular mass and amino acid sequence 
information. These two technologies have been proved to be 
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complementary for a comprehensive comparative proteomic 
analysis (15, 16). ICAT analysis showed a clear bias for pro- 
teins with high molecular mass, whereas the 2D electrophoresis 
or DIGE method could separate proteins in certain low molec- 
ular mass ranges and also identified cysteine-free proteins that 
were transparent to the ICAT analysis. Moreover ICAT analysis 
quantifies the sum of the protein species of one gene product, 
whereas the 2D electrophoresis or DIGE method quantifies at 
the level of resolved protein species, including post-transla- 
tionally modified and processed polypeptides. 

In the present work, systematic analyses of the proteome of 
SARS-CoV-infected cells were performed using 2D-DIGE fol- 
lowed by electrospray mass spectrometry identification and 
ICAT technology coupled with 2D-LC-MS/MS. Using these 
two complementary methods, the differentially expressed 
proteome profiles between Vero E6 cells with and without 
virus infection were created, and a total of 355 proteins or 
protein spots were identified and quantitated. Further analysis 
of these data provides the clues for understanding the infec- 
tion and pathogenesis of the SARS-CoV and the virus-host 
interactions. 


EXPERIMENTAL PROCEDURES 


Materials— Analytical reagent grade chemicals were used through- 
out unless otherwise stated. Water was prepared using a Milli-Q water 
purification system (Millipore, Bedford, MA). Chemicals used for gel 
electrophoresis were purchased from Bio-Rad. Cy2, Cy3, Cy5, Pre- 
cast IPG dry strips, pH 3-10 non-linear, and ECL PLUS Western 
blotting detection regents were purchased from Amersham Bio- 
sciences. Formic acid and guanidine hydrochloride were obtained 
from Sigma. HPLC grade ACN was from Fisher. Sequencing grade 
trypsin was obtained from Promega (Southampton, UK). 

Cell Culture, Virus Infection, and Sample Preparation—African 
green monkey kidney cells (Vero E6, ATCC) were maintained in Dul- 
becco’s modified Eagle’s medium (DMEM; Invitrogen) supplemented 
with 10% fetal bovine serum (FBS; Invitrogen) at 37 °C in 5% COs. 

For virus infection, Vero E6 cells were treated with the DMEM (2% 
FBS) containing SARS-CoV virions (BJ-01 isolate, provided by Acad- 
emy of Military Medical Sciences) for 1 h. The virus medium was 
removed after the infection, and the infected cells were cultured in 
DMEM with 2% FBS at 37 °C in 5% COs. 24 h postinfection, the cells 
were washed with PBS twice and then lysed with a lysis buffer 
containing 40 mm Tris base, 60 mm DTT, 8 M urea, 4% CHAPS. The 
mixture was centrifuged in a microcentrifuge for 5 min, and the 
supernatant was collected as the infected cellular sample. All the 
experiments using the virus were carried on in a biosafety Level 3 
laboratory. 

DIGE—For DIGE analysis, cellular samples were precipitated over- 
night with 5 volumes of 50:50:0.1 volumes of ethanol:acetone:acetic 
acid at —20 °C and resolubilized in the lysis buffer (7 m urea, 2 M 
thiourea, 30 mm Tris-Cl, 4% CHAPS). The cell lysates (E6 or E6-V) 
were labeled with N-hydroxysuccinimidyl ester derivatives of the cya- 
nine dyes Cy2, Cy3, and Cy5 (Amersham Biosciences) according to 
the manufacturer’s guidelines. Typically 50 wg of lysate were mini- 
mally labeled with 400 pmol of Cy2, Cy3, or Cy5, respectively. Differ- 
entially labeled samples were mixed (Cy2-labeled standard, Cy3- 
labeled E6, and Cy5-labeled E6-V for one gel; Cy2-labeled standard, 
Cy5-labeled E6, and Cy3-labeled E6-V for another gel). An equal 
volume of 2x 2D sample buffer (7 m urea, 2 m thiourea, 4% CHAPS, 
2% DTT, and 2% IPG buffer, pH 3-10 non-linear) was added to 


provide 65 mm DTT to reduce the labeled samples for 15 min at room 
temperature. Then the samples were subjected to 2D-PAGE in the 
dark. The DIGE gels were scanned using Typhoon Variable Model 
Imagels 9400 (Amersham Biosciences). DeCyder (Amersham Bio- 
sciences) software was used for image analysis. 

In-gel Trypsin Digestion—Changes observed in 2D-DIGE images 
were aligned with Phastgel Blue R-stained protein patterns in prepar- 
ative gels. Spots of interest were cut manually. Gel pieces were 
digested as described by Yu et al. (17). 

RP-HPLC-ESI-MS/MS Analysis— RP-HPLC was performed using a 
surveyor LC system (Thermo Finnigan, San Jose, CA) on a Cg col- 
umn (RP, 180 wm x 150 mm, BioBasic® C,,, 5 wm, Thermo Hypersil- 
Keystone). The pump flow rate was split 1:120 to achieve a column flow 
rate of 1.5 ul/min. Mobile phase A was 0.1% formic acid in water, and 
mobile phase B was 0.1% formic acid in acetonitrile. The tryptic peptide 
mixtures were eluted using a gradient of 2-98% B over 60 min. 

The mass spectral data shown here were acquired on an LTQ linear 
ion trap mass spectrometer (Thermo Finnigan) equipped with an 
electrospray interface operated in positive ion mode. The temperature 
of the heated capillary was set at 170 °C. A voltage of 3.4 kV applied 
to the ESI needle resulted in a distinct signal. Normalized collision 
energy was 35.0. The mass spectrometer was set so that one full MS 
scan was followed by three MS/MS scans on the three most intense 
ions from the MS spectrum with the following Dynamic Exclusion™ 
settings: repeat count, 2; repeat duration, 0.5 min; exclusion duration, 
2.0 min. 

ICAT Analysis—|CAT analysis was performed using the Cleavable 
ICAT™ reagent kit (Applied Biosystems, Foster City, CA) according to 
the manufacturer’s guidelines with some modifications. For ICAT 
analysis, the cellular samples were precipitated and resolubilized in 
denaturing buffer (6 m guanidine hydrochloride, 100 mm Tris-Cl, pH 
8.3). 100 ug of the E6 or E6-V protein sample in 80 yl of denaturing 
buffer were reduced at 37 °C for 2 h with 5 mm tributylphosphine 
(Bio-Rad) and alkylated at 37 °C for 2 h in the dark with ICAT-light and 
ICAT-heavy reagent, respectively. The samples were digested with 
trypsin at 37 °C for 20 h. Then the ICAT-labeled peptides were puri- 
fied using the kit of ICAT™ Avidin Buffer Pack and Avidin Affinity 
Cartridge (Applied Biosystems) according to the manufacturer’s 
guidelines. 

2D LC-MS/MS — Orthogonal 2D LC-MS/MS was performed using a 
ProteomeX work station (Thermo Finnigan). The system was fitted 
with a strong cation exchange column (320-m inner diameter x 100 
mm, DEV SCX, Thermo Hypersil-Keystone) and two C,g RP columns 
(180 xm x 100 mm, BioBasic C,,, 5 wm, Thermo Hypersil-Keystone). 
The salt steps used were 0, 25, 50, 75, 100, 150, 200, 400, and 800 
mm NH,Cl synchronized with nine 140-min RP gradients. RP solvents 
were 0.1% formic acid in either water (A) or acetonitrile (B). The 
setting of the LCQ Deca Xplus ion trap mass spectrometer is as 
follows. One full MS scan was followed by three MS/MS scans on the 
three most intense ions from the MS spectrum according to the 
following Dynamic Exclusion settings: repeat count, 1; repeat dura- 
tion, 0.5 min; exclusion duration, 3.0 min. 

Data Base Searching—The acquired MS/MS spectra were auto- 
matically searched against the nonredundant human protein data 
base (NCBI (www.ncbi.nlm.nih.gov), December 4, 2003 release) using 
the TurboSEQUEST program in the BioWorks™ 3.1 software suite. 
For protein spot identification, an accepted SEQUEST result had to 
have a ACn score of at least 0.1 (regardless of charge state). Peptides 
with a +1 charge state were accepted if they were fully tryptic and 
had a cross correlation (Xcorr) of at least 1.8. Peptides with a +2 
charge state were accepted if they had an Xcorr >2.5. Peptides with 
a +3 charge state were accepted if they had an Xcorr >3.7. For ICAT 
analysis, protein identification and quantification were achieved by 
using SEQUEST and EXPRESS software tools. Peptides with a +1 
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Fig. 1. Cy dye image of the SARS- 
CoV-infected Vero E6 (red, Cy3-la- 
beled) cells and uninfected Vero E6 
cells (green, Cy5-labeled). The image 
shown is of a 100-ug sample (60 pg 
each of Cy3- and Cy5-labeled lysates) 
run on a pH 3-10 non-linear gradient IPG 
strip and 12.5% polyacrylamide gel. The 
red or green spots indicate the differen- 
tially expressed proteins. The spots cir- 
cled and marked with numbers have 
been identified and listed in Supplemen- 
tal Table |. 
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TABLE | 
The number of differentially expressed proteins in SARS-CoV-infected or uninfected Vero E6 cells identified by DIGE and ICAT analysis 


Ratio of E6:E6-V =<0.50 <0.67 = 0.83 0.83-1.20 =1.20 =1.50 =2.00 Total 
DIGE (unique proteins/protein spots) 11/14 17/21 23/28 0/0 27/35 13/15 3/3 48/63 
ICAT (unique proteins) 37 48 69 48 205 119 51 322 
Total (unique proteins) 46 60 86 48 221 126 52 355 


charge state were accepted if they were fully tryptic and had an Xcorr 
of at least 1.5. Peptides with a +2 charge state were accepted if they 
had an Xcorr >2.0. Peptides with a +3 charge state were accepted if 
they had an Xcorr >2.5. Then the confirmation of protein identification 
and quantification of the peptides was further analyzed manually as 
described by Han et al. (14). Protein abundance ratios larger than 
+1.5 or smaller than —1.5 (or 0.67) were set as a threshold indicating 
significant changes (18). 

Western Blotting —SDS-PAGE-separated proteins were transferred 
onto a nitrocellulose membrane (Protran, Schleicher & Schuell) on a 
Mini Trans-Blot Cell (Bio-Rad). The transferred membrane was 
blocked with 1 NET-Gelatin (150 mm NaCl, 5 mm EDTA, pH 8.8, 50 
mm Tris-HCl, 0.05% Triton X-100, 8.3% gelatin) at room temperature 
for 1 h and then was incubated with corresponding primary antibodies 
(anti-HSP90 (rabbit), Santa Cruz Biotechnology, 1:5,000; anti-a-tubu- 
lin (mouse), Sigma, 1:600,000; anti-stathmin (rabbit), Calbiochem- 
Novabiochem, 1:5,000; anti-vimentin (mouse), Sigma, 1:1,000) for 1h 
at room temperature, and then the membranes were incubated with 
the corresponding secondary antibodies. After being detected with 
ECL PLUS (Amersham Biosciences) according to the manufacturer’s 


instructions, the membrane was scanned using a Typhoon Scanner 
9400 (Amersham Biosciences). 

Bioinformatic Annotation—The theoretical pl and molecular mass 
values of proteins were defined by the program pepstats (www. 
hgmp.mrc.ac.uk/Software/EMBOSS). The protein function and sub- 
cellular location annotation was from the Swiss-Prot and TrEMBL 
protein data base (us.expasy.org/sprot/). 


RESULTS 


2D-DIGE Analysis of the SARS-CoV-infected and Unin- 
fected Vero E6 cells—2D-DIGE as a qualitative and quantita- 
tive proteomic approach was performed to determine the 
differential proteomes of the SARS-CoV-infected and unin- 
fected Vero E6 cells. The lysate from infected cells was la- 
beled with Cy5, whereas the lysate from uninfected cells was 
labeled with Cy3. After electrophoresis and imaging, the Cy3 
and Cy5 images were false colored in red and green, respec- 
tively, and two images were overlapped (Fig. 1). The biological 
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TABLE II 
List of protein spots with at least 2.0-fold quantitative alteration in SARS-CoV-infected or uninfected Vero E6 cells based on DIGE analysis 


VAMP, vesicle-associated membrane protein; PP2A, protein phosphatase 2A. 


; : : ‘ Ratio 
F 5 dt Subcellular Theoretical Experimental Peptide Unique Sequence ? 
DIGE Sporn: RYOTE GeSchplion Ging. location MM (Da)/pl MM (kDa)/pl hits peptides coverage (SEEM) 
% 
Actin network protein 
42 Calponin, acidic gi|6225157| Cytoskeletal 36,413.61/5.69  36.5/5.77 12 8 27.05 1:0.38 
isoform 
Chromosome-associated protein 
20 Chromatin gil1172846] Nuclear 47,655.74/4.74 = 55.2/4.71 4 4 13.18 1:0.34 
assembly factor 1 
subunit C 
Cytoskeletal protein 
21 a-Tubulin 2 gil20455316| Cytoskeletal 50,151.71/4.94  55.4/5.20 15 12 39.69 1:2.00 
Metabolic enzyme protein 
18 Pyruvate kinase, gil20178296| No annotation 57,913.83/7.96  34.3/5.89 25 10 19.21 1:0.37 
M1 isozyme 
19 Pyruvate kinase, gil20178296| No annotation 57,913.83/7.96  34.6/7.63 18 10 20.72 1:0.38 
M1 isozyme 
38 Galactokinase gil1730187| No annotation 42,272.43/6.04 38.9/5.76 11 6 13.78 1:0.32 
43 Transaldolase gi|6648092| Cytoplasmic 37,540.16/6.36  36.0/6.34 5 4 13.35 1:0.40 
44 Pyruvate gi|129070| Mitochondrial 39,219.39/6.20 31.4/5.43 18 9 29.25 1:0.49 
dehydrogenase 
E1 component B 
subunit 
Heterogeneous nuclear ribonucleoprotein 
32 hnRNP F gil1710628] Nuclear 45,671.98/5.38 47.6/5.22 16 10 33.01 1:2.28 
Intracellular trafficking protein 
12 Transmembrane gil19920317| Membrane 66,022.66/5.63 62.0/4.76 7 7 15.61 1:0.27 
protein (63 kDa) 
13 Transmembrane gil19920317| Membrane 66,022.66/5.63  61.7/4.78 5 4 10.30 1:0.34 
protein (63 kDa) 
52 VAMP-associated gil24638339| Membrane 27,228.36/6.85 27.3/7.49 5 5 21.81 1:0.49 
protein B/C 
Protein processing protein 
29 Mitochondrial gil29840827| Membrane 54,366.18/6.38 48.4/5.84 9 5 11.86 1:0.40 
processing 
peptidase B 
subunit 
Signaling protein 
14 PP2A subunit A, gil231443| No annotation  65,223.6/4.96 59.8/4.89 7 6 11.04 1:0.31 
PR65-a isoform 
47 14-3-3 protein 7 gi|112690| Cytoplasmic 27,764.29/4.68  27.5/6.52 3 2 8.98 1:2.51 
48 Annexin IV gi|1703319| No annotation 35,882.74/5.84  28.0/5.47 6 20.06 1:0.27 
55 Growth factor gi|121603| No annotation 25,206.38/5.89 25.2/5.90 11 8 31.80 1:0.28 
receptor-bound 
protein 2 


variation analysis mode of DeCyder was used for comparing 
the overlay image. The standardization was done by the com- 
parison of normalized Cy3 and Cy5 protein spot volumes with 
the corresponding Cy2 standard spot volumes within each 
gel. Protein spots with an average ratio value greater than 
1.2-fold and at test p value <0.05 were selected for mass 
spectrometric identification. A total of 63 proteins belonging 
to 48 unique gene products were identified with ESI-LC- 
MS/MS (Supplemental Table |). Among those proteins, 17 
proteins (21 protein spots) were down-regulated, whereas 13 
proteins (15 protein spots) were up-regulated by at least 
1.5-fold in SARS-CoV-infected cells (Table I). The identified 
protein spots with at least a 2.0-fold quantitative alteration in 
SARS-CoV-infected Vero E6 cells are listed in Table Il. The 


differentially expressed proteins are involved in various func- 
tions including cytoskeleton, actin-associated network, met- 
abolic enzymes, and signal transduction (Table II and Supple- 
mental Table 1). 

ICAT Analysis of the SARS-CoV-infected and Uninfected 
Vero E6 Cells—Recently a new proteomic approach, ICAT, 
has been developed to quantitatively analyze the protein dif- 
ferential expression, which is faster than the gel method. In 
the present work, the ICAT method coupled to 2D-LC-MS/MS 
was used for quantitative comparison of differential proteome 
profiles between SARS-CoV-infected cells and uninfected 
cells. Cysteines were labeled with light ('*C) and heavy ('C) 
ICAT reagent, and labeled peptides were affinity-purified us- 
ing an avidin column. The resultant peptides were first frac- 


Molecular & Cellular Proteomics 4.7 905 


ASBMB 


ow 
a 
~~ 


pe, 


MOLECULAR & CELLULAR PROTEOMICS 


MCP 


Quantitative Analysis of SARS Coronavirus-infected Cells 


TABLE III 
Selected list of proteins in SARS-CoV-infected or uninfected Vero E6 cells based on ICAT analysis 


L, light; H, heavy; MuSK, muscle-specific tyrosine kinase; NSF, N-ethylmaleimide-sensitive factor. 


: ‘at Subcellular Peptide Unique Ratio, ICAT Peptide sequence 
Protein description GI no. location hits peptides (L:H; E6:E6-V) (charge/Xcorr/ACn)? 
Actin network protein 
Anillin, actin-binding protein gil31657094| Cytoskeletal 1 1 1:0.25 R.FGERC*QEHSKESPAR:S (2+/2.05/0.20) 
a-Actinin 3 (F-actin cross- gi|728751| Cytoskeletal 1 1 1:0.38 K.C*QLEINFNTLQTK.L (2+/3.16/0.36) 
linking protein) 
Destrin (actin-depolymerizing gi|5802966| Cytoskeletal 4 1 2.77 K.HEC*QANGPEDLNR.A (3+/3.34/0.37) 
factor) 
Channel or transporter protein 
Potassium channel-modulatory gi/31543383] No annotation 1 1 1:2.59 R.HEGVSC*DAC*LK.G5 (2+/2.90/0.43) 
factor 1 
Potassium channel subfamily K gi/24636282| Membrane 1 1 1:2.97 K.VFC*MFYALLGIPLTLVTFQSLGER.L 
member 15 (3+/2.52/0.11) 
Sodium bicarbonate gi]15042959] No annotation 1 1 1:8.84 R.FFGGLC*LDIK.R (2+/2.19/0.20) 
transporter 4 isoform a 
Chromosome-associated protein 
DEK protein gi|544150| Nuclear 1 1 1:12.45 K.C*PEILSDESSSDEDEKK.N (3+/2.90/0.20) 
High mobility group protein 4 gil20138144] Nuclear 1 1 1:2.20 K.RPPSGFFLFC*SEFRPK.| (3+/3.70/0.34) 
Condensin subunit 2 gi/30172801| Cytoplasmic 1 1 1:2.21 K.TAASFDEC*STAGVFLSTLHC*QDYR.S 
and nuclear (3+/4.73/0.54) 
Heat shock protein 
HSP 90-a gil123678| Cytoplasmic 3 1 1:2.08 K.HGLEVIYMIELIDKYC*VQQLK.E 
(3+/3.58/0.11) 
Stress-induced gi|400042| No annotation 2 1 1:2.18 K.DC*EEC*IQLEPTFIK.G (2+/2.54/0.25) 
phosphoprotein 1 
Heterogeneous nuclear 
ribonucleoprotein® 
hnRNP DO gil13124489| Nuclear 4 1 1151 R.GFC*FITFK.E (2/2.27/0.41) 
hnRNP U gil6226894| Nuclear 4 1 1:1.77 K.TC*NC*ETEDYGEK.F (2+/2.68/0.37) 
hnRNP core protein A1 gil133254| Nuclear 6 1 1:1.90 K.YHTVNGHNC*EVR.S (2+/3.58/0.24) 
hnRNP | gil131528| Nuclear 2 1 1:1.95 K.LSLDGQNIYNAC*C*TLR.1 (2+/3.37/0.41) 
hnRNP K isoform a gil14165437| Cytoplasmic 4 1 1:2.00 K.LFQEC*C*PHSTDR.V (2+/3.55/0.45) 
and nuclear 
hnRNP AO gi/8134660 Nuclear 1 1 1:2.04 R.GHFEAFGTLTDC*VVVVNPQTK.R 
(3+/2.98/0.43) 
hnRNP D-like gil4885423 No annotation 8 1 1:2.18 R.RGFC*FITYTDEEPVKK.L (3+/4.30/0.43) 
hnRNP E2 gil6707736 Nuclear 1 1 1:2.27 R.INISEGNC*PER.| (2+/2.84/0.29) 
Intracellular trafficking protein 
Translocation protein SEC63 gi|18203500| Endoplasmic 1 1 1:0.50 R.APTLASLENCMKLSQMAV- 
homolog reticulum QGLQQFK.S (3+/2.74/0.29) 
Synaptophysin-like protein gi/5803185| Vesicular 1 1 1:2.28 K.GQTEIQVNC*PPAVTENKTVTATFGYPFR.L 
isoform a; pantophysin (3+/2.71/0.21) 
a-Soluble NSF attachment gil6226705 Membrane 1 1 1:9.51 K.C*LLKVAGYAALLEQYQK.A 
protein (3+/2.66/0.13) 
Metabolic enzyme 
Vacuolar ATP synthase gil22096378| Endoplasmic 1 1 1:0.36 R.FCPFYKTVGMLSNMIAFYDMARR.A 
catalytic subunit A reticulum (3+/3.26/0.21) 
Pyruvate dehydrogenase E1 gi|129063| Mitochondrial 1 1 1:0.36 K.LPCIFICENNR.Y (2+/2.13/0.23) 
component a subunit 
FK506-binding protein 4 gi|399866| Cytoplasmic 2 1 1:0.41 K.VGEVCHITCKPEYAYGSAGSPPK.| 
and nuclear (3+/2.86/0.23) 
Alanine-glyoxylate gil17432913] Mitochondrial 1 1 1:0.45 K.CLQHFNTFGGNPMACAIGSAVLEVIK.E 
aminotransferase 2 (3+/2.77/0.16) 
GMP reductase 2 gil25008511| No annotation 1 1 1:2.10 K.GHIISDGGC*SC*PGDVAK.A 
(3+/2.63/0.16) 
Succinyl-CoA synthetase, a gil20141765| Mitochondrial 6 1 1:2.16 R.LIGPNC*PGVINPGEC’*K.| (2+/3.78/0.52) 
chain 
Thioredoxin-like protein p46 gil29839560| Endoplasmic 19 2 1:2.44+0.45 K.VDC*TAHSDVC*SAQGVR.G (2+/4.18/ 
reticulum 0.50); K.VDC*TQHYELC*SGNQVR.G 
(3+/4.27/0.35) 
UDP-glucose 6-dehydrogenase gi|6175086| No annotation 1 1 1:3.56 R.AVQALC*AVYEHWVPR.E (3+/2.60/0.30) 
Protein degradation 
Ovochymase gil34419641| No annotation 1 1 1:2.03 R.YLLDYRGRLEC*SWVLR.V (3+/2.62/0.17) 
Proteasome subunit P50 gil20532406| Cytoplasmic 3 1 1:2.22 R.C*TDDFNGAQC’*K.A (2+/3.16/0.45) 
and nuclear 
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TABLE III— continued 


Protein description GI no. ae rae peur ear eeee y Peptide sequence (charge/Xcorr/ACn)* 
Ribosomal protein? 
Ribosomal protein S4 gil4506725 Ribosomal 11 2 1:1.52 + 0.33 K.LREC*LPLIIFLR.N (8+/3.82/0.32) 
60 S ribosomal protein L10a gil15431288| Ribosomal 6 1 1:1.58 K.FSVC*VLGDQQHC*DEAK.A 
(2+/2.97/0.37) 
40 S ribosomal protein S15a gil14165469| Ribosomal 1 1 1:1.61 K.C*GVISPR.F (2+/2.10/0.16) 
40 S ribosomal protein S21 gi|4506699 Ribosomal 2 1 1:1.65 K.TYAIC*GAIR.R (2 +/2.27/0.25) 
40 S ribosomal protein S26 gil15011936| Ribosomal 6 1 1:1.66 K.LHYCVSCAIHSK.V (8+/3.10/0.53) 
40 S ribosomal protein S12 gil14277700| Ribosomal 9 2 1:1.66 + 0.12 K.LVEALC*AEHQINLIK.V (2+/4.17/0.19); 
R.KVVGC*SC*VVVK.D (2+/2.72/0.21) 
40 S ribosomal protein S11 gil|4506681 Ribosomal 13 1 1:2.03 R.DVQIGDIVTVGEC*RPLSK.T 
(2+/4.65/0.53) 
40 S ribosomal protein SA gi|125969| Ribosomal 4 1 1:2.06 R.ADHQPLTEASYVNLPTIALC*NTDSPLR.Y 
(3+/5.15/0.49) 
40 S ribosomal protein S27a gil4506713 Ribosomal 7 2 1:2.18 + 0.23 K.C*C*LTYC*FNKPEDK.- (2+/3.35/0.40); 
R.EC*PSDEC*GAGVFMASHFDR.H 
(3+/4.17/0.45) 
Signaling protein 
Cell division protein kinase 6 gil266423| No annotation 1 1 1:0.17 R.ADQQYECVAEIGEGAYGK.V 
(3+/2.51/0.13) 
Ectodysplasin A receptor- gil21362527| Cytoplasmic 1 1 1:0.40 K.ENCTCSSCLLR.A (2+/2.25/0.15) 
associated adapter protein 
Transducin B chain 5 gi|38258891| No annotation 1 1 1:0.45 R._YYPSGDAFASGSDDATCRLYDLR.A 
(3+/2.51/0.22) 
Phosphoinositide-3-kinase, gil23943912| No annotation 1 1 1:0.49 K.PVIPVLSSTILPSTY QIRITTO*K.T 
regulatory subunit 4, p150 (3+/2.55/0.12) 
Galectin-1 gil126155| No annotation 23 2 1:2.07 + 0.35 K.DSNNLC*LHFNPR.F (2+/3.34/0.32); 
R.FNAHGDANTIVO*NSK.D (2+/4.70/0.59) 
Receptor tyrosine kinase gil5031927| No annotation 1 1 1:2.18 R.EYC*LAVKELFC*AKEWLVMEEK.T 
MuSK (3+/2.82/0.15) 
Fibroblast growth factor-9 gi|544290| Secreted 1 1 1:2.42 R.QLYCRTGFHLEIFPNGTIQGTR.K 
(3+/2.78/0.22) 
Apoptotic chromatin gil7662238| Nuclear 1 1 1:2.47 K.FLC*ADYAEQDELDYHR.G (3+/3.36/0.42) 
condensation inducer in the 
nucleus 
Testin gil17380320| No annotation 2 1 1:2.52 K.NHAVVC*QGC*HNAIDPEVQR.V 
(3+/3.03/0.37) 
Tenascin X precursor gil9087217 Secreted 1 1 1:3.68 R.VRGEESEVTVGGLEPGC*K.Y 
(3+/3.58/0.1 1) 
Transcription or replication factor 
DNA ligase | gi|118773| Nuclear 1 1 1:0.38 K.GLFVACRHSEARFIAR.S (3+/2.65/0.22) 
Cleavage- and gil5901928 Nuclear 1 1 1:2.01 K.RELHGQNPVVTPC*NK.Q (3+/2.69/0.40) 
polyadenylation-specific 
factor 6, 68-kDa subunit 
DNA helicase homolog gil5523990 Nuclear 1 1 1:4.31 R.QQLPLQLAWAMSIHKSQGMTLDC*VEISL- 
GR.V (3+/3.08/0.12) 
Translation factor? 
elF-1 gil1174483 No annotation 8 2 1:1.56 + 0.08 K.FAC*NGTVIEHPEYGEVIQLQGDQR.K 
(3+/4.16/0.38); KINICQFLVEIGLAK.D 
(2+/2.20/0.15) 
elF-5 gil27735202| No annotation 1 1 1:1.60 K.FVLCPECENPETDLHVIPK.K (3+/2.93/ 
0.27) 
elF-5A2 protein gil9966867| No annotation 1 1 1:1.76 K.KYEDIC*PSTHNMDVPNIK.R (3+/3.42/ 
0.41) 
EF-Tu, mitochondrial precursor gil1706611| | Mitochondrial 4 1 1:2.17 R.HYAHTDC*PGHADYVK.N (3+/3.96/0.49) 
Other function protein 
WD repeat domain 17 isoform — gi31317311| No annotation 1 1 1:0.42 R.IWDYTQDACINILNGHTAPVR.G 
1 (3+/2.52/0.15) 
Ran-binding protein 2 gil1709217| Nuclear 3 1 1:2.03 K.C*IAC*QNPGK.Q (2+/2.86/0.43) 
Hyaluronan-binding protein 4 gil24307947| No annotation 1 1 1:2.06 R.YGGNDKIAVRTEDNMGGC*GVR.T 
(3+/2.60/0.17) 
Metallothionein II gi|127397| No annotation 1 1 1:2.51 K.CAQGCICK.G (2+/2.20/0.15) 


? Select list of the peptides with the best assignment scores for protein identification and quantification by ICAT analysis. C/C*, ICAT 
light/heavy reagent-labeled cysteine-including peptide, respectively. 

© For heterogeneous nuclear ribonucleoproteins, ribosomal proteins, and translation factors, proteins with at least 1.5-fold alteration are 
listed. The MS and MS/MS spectra for identification and quantitation of these proteins are shown in Supplemental Figs. 4-24. 
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tionated with nine salt steps, and each step was further sep- 
arated using a reverse phase column and analyzed with ESI- 
MS/MS. The relative quantitation of proteins in infected and 
uninfected cells was determined from the relative abundance 
ratio of labeled peptides. By this approach, the relative abun- 
dance of 322 proteins was measured (Supplemental Table 1). 
Individual peptide sequences and their probability scores for 
identification and quantitation of those proteins are listed in 
Supplemental Table Il. Among these cellular proteins identi- 
fied from the infected cells, 167 proteins had at least 1.5-fold 
quantitative alterations with 119 proteins up-regulated and 48 
proteins down-regulated (Table l). Some selected proteins 
with quantitative alterations are listed in Table Ill. In addition, 
two SARS proteins, M and S, were observed to be signifi- 
cantly increased in the infected cells (data not shown). 

Comparison and Validation of the Quantitation of Differential 
Proteins—We used 2D-DIGE and ICAT methods to investi- 
gate the differentially expressed proteins in Vero E6 cells 
infected with SARS-CoV. These two methods are both more 
sensitive and accurate than traditional methods such as silver 
staining (data not shown). We further used Western blot assay 
to confirm the differential expression of the proteomes iden- 
tified by those two methods. The Western blotting results 
showed that the ratios of four representative proteins, a-tu- 
bulin, HSP90, OP18 stathmin, and vimentin, between the 
infected and uninfected cells were in agreement with those 
obtained from 2D-DIGE or ICAT approaches (Fig. 2 and Sup- 
plemental Figs. 1-3). 

It was observed that the overall expression alterations were 
similar with both quantitative proteomic approaches, although 
only 15 differential proteins were identified by both 2D-DIGE 
and ICAT approaches (Table IV). In addition, the DIGE method 
detected more isoforms of proteins than the ICAT approach 
did (Table IV). 

Subcellular Location and Function Classification of the Dif- 
ferential Proteins—\n the present work, a total of 355 unique 
gene products of SARS-CoV-infected and uninfected cells 
were identified and quantitated with either the 2D-DIGE or 
ICAT method of which 186 proteins had at least 1.5-fold 
quantitative alterations (Table | and Supplemental Table |). 
Among those 186 differentially expressed proteins, 60 pro- 
teins were down-regulated in infected cells, whereas 126 
proteins were up-regulated. These proteins were further clas- 
sified according to their subcellular locations (Fig. 3). Interest- 
ingly, the up-regulated proteins in infected cells were mainly 
located in the nuclei (about 25%, see Fig. 3A), whereas down- 
regulated proteins distributed within the cells (Fig. 3B). 

Those 186 proteins with at least 1.5-fold differential expres- 
sion were further classified according to their functions (Sup- 
plemental Table |). Fig. 4 presents the functional categories of 
these proteins. In these identified SARS-CoV-infected cellular 
proteins, the enzymes, signal proteins, ribosomal proteins, 
and heterogeneous nuclear ribonucleoproteins (hnRNPs) con- 
sisted of about 45% of the total up-regulated proteins (Fig. 


E6 E6-V E6 E6-V 


I DIGE 1:2.00 ICAT tt. DIGE 1:1.51 BCAT 1:1.) 


m -- -— — “— 


A a- tubulin B HSP90 
E6 E6-V E6 E6-V 
I DIGE 1:1.62 CAT BLA? 


C Stathmin (OP18) D Vimentin 


Fic. 2. Protein quantitative confirmation with Western blotting. 
A, a-tubulin; B, HSP90; C, OP18 stathmin; D, vimentin. Row / indi- 
cates the DIGE or ICAT analysis ratio. Row // shows the three- 
dimensional fluorescence intensity profiles of the individual spots. 
Row III shows the Western blot results. The MS and MS/MS spectra 
of ICAT analysis for proteins a-tubulin, HSP90, and vimentin are 
shown in Supplemental Figs. 1-3, respectively. 


4A). On the other hand, more than half of the down-regulated 
proteins of the infected cells were involved in the enzymatic 
reactions, signal transduction, immune responses, and actin 
networks (Fig. 4B). 


DISCUSSION 


From the literature, very few studies have been performed 
to analyze the interaction between coronavirus and host cells 
with proteomic approaches. In the present work, we used two 
quantitative proteomic assays, DIGE and ICAT, to determine 
the differentially expressed protein profiles of SARS-CoV- 
infected and uninfected cells. It was noted that only a few 
proteins were identified by both DIGE and ICAT approaches, 
suggesting that these two methods are complementary to 
each other. In addition, the main difference of those two 
assays is that ICAT only determines the overall expression 
level, whereas DIGE can detect and quantitate protein iso- 
forms in a gel. Although the DIGE method can detect protein 
isoforms possibly caused by post-translational modifications, 
it has limitations in identification of proteins with very high or 
low molecular weights, extreme acidic/alkaline proteins, and 
low abundance proteins. In our present work, 45 proteins that 
have a molecular mass greater than 100 kDa were identified 
among which only two proteins were detected by the DIGE 
approach. Moreover all 22 identified proteins with pl > 10 
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TABLE IV 
List of the proteins or protein spots identified and quantitated in SARS-CoV-infected or uninfected Vero E6 cells by both ICAT and DIGE analysis 


Ratio 
F 4 at Subcellular Theoretical DIGE Experimental Peptide Unique Sequence 
Protein description Gling. location MM (Da)/p!|__ spot no. MM (kDa)/pl hits peptides paver ICAT (L:H; DIGE 
E6:E6-V) —_ (E6:E6-V) 
% 
Actin network protein 
Cofilin, non-muscle gil116848] Cytoplasmic 18,502.47/8.22 60 19.0/8.30 14 2 15.66 1:1.59 + 0.08 
isoform and nuclear 61 19.1/6.82 1 1 8.43 1:1.81 
120 35 56.12 1:0:85 
Chromosome- 
associated protein 
Chromatin assembly —_gil1172846| Nuclear 47,655.74/4.74 1 1 3.76 1:1.63 
factor 1 subunit C 20 55.2/4.71 4 4 13.18 1:0.34 
Cytoskeletal protein 
a-Tubulin 2 gi|20455316| Cytoskeletal 49,959.62/4.98 2 2 8.22 1:1.90 + 0.10 
21 55.4/5.20 15 12 39.69 1:2.00 
Heat shock protein 
HSP 90-8 gil17865718| Cytoplasmic 83,264.38/4.97 1 1 2.62 1:1.23 
6 82.5/5.31 54 24 34.53 
5 82.9/5.27 20 14 26.80 
Heat shock cognate gi|123648| No annotation 70,897.96/5.37 7 1 3.41 1:1.34 
71-kDa protein 9 70.8/5.78 13 10 19.04 1:0.52 
11 70.4/4.18 15 11 21.05 1:0.54 
10 71.0/4.00 33 19 30.65 1:0.72 
Hsc70-interacting gi|6686278| Cytoplasmic 41,331.83/5.18 3 1 4.61 1:1.41 
protein 31 48.0/4.99 13 v4 15.99 1:1.73 
Stress-induced gi|400042| No annotation 62,639.39/6.40 2 1 2.58 1:2.18 
phosphoprotein 1 15 63.4/7.35 81 30 46.78 1:1.22 
Heterogeneous nuclear 
ribonucleoprotein 
hnRNP H gil1710632| Nuclear 49,229.59/5.89 1 1 4.23 1:1.57 
28 48.9/6.04 29 12 33.41 1:1.68 
27 51.0/5.83 19 11 34.97 1:1.60 
26 51.0/5.83 49 12 30.96 1:1.46 
Metabolic enzyme 
protein 
Nucleoside gil127983| Cytoplasmic, 17,298.03/8.52 1 1 5.92 1:1.05 
once ne kinase and nuclear 62 13.3/7.92 2 1 5.92 1:0.57 
a enolase gil119339| Cytoplasmic 47,168.91/7.01 19 3 11.98 1:1.48 + 0.08 
33 47.2/7.00 136 25 50.00 0 
34 47.4/7.67 188 20 46.08 
Pyruvate kinase, M1 gi/20178296| No annotation 57,936.87/7.96 17 5 11.30 1:1.75 + 0.44 
isozyme 22 57.5/8.21 120 35 56.12 1:1.42 
17 58.0/8.00 9 v4 12.43 1:1.22 
18 34.3/5.89 25 10 19.21 1:0.37 
19 34.6/7.63 18 10 20.72 1:0.38 
Thioredoxin gi|135773| No annotation 11,737.48/4.82 5 1 8.57 1:1.77 
63 12.0/4.80 7 5 27.62 1:1.28 
UDP-glucose 6- gil6175086| No annotation 55,024.12/6.73 1 1 3.04 1:3.56 
dehydrogenase 23 56.6/8.08 14 13 29.35 1:1.27 
Signaling protein 
14-3-3 protein 7 gil112690| Cytoplasmic 27,764.29/4.68 3 1 4.49 1:1.57 
47 27.5/6.52 3 2 8.98 
46 27.8/4.60 2 1 5:71 9 
Annexin IV gi|1703319| No annotation 35,882.74/5.84 5 2 7.52  1:1.62 + 0.19 
48 28.0/5.47 9 6 20.06 1:0.27 


were contributed by the ICAT method (Supplemental Table |). 
Concerning the capacity of the protein identification, the DIGE 
approach identified 48 unique proteins, whereas the ICAT 
approach obtained 322 proteins; the ICAT approach espe- 
cially identified more signal proteins, which usually are low 
abundance proteins in a cell. 

Based on the identified proteins in the present work, we can 
gain an overall insight into the altered protein expression of the 


host cell responding to SARS-CoV infection. Among the differ- 
entially expressed host proteins, many of them participated in 
viral RNA replication and translation (Tables II and III and Sup- 
plemental Table 1). It is known that positive-strand RNA viruses 
recruit normal components of host cellular RNA-processing or 
translation machineries for the viral RNA synthesis and protein 
synthesis (19, 20). Our data here indicate that SARS-CoV virus 
uses a similar strategy when infecting the host cells. 
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Up-regulated proteins in SARS-CoV infected cells with at least 
1.5 fold quantitative alterations 


Fic. 3. Subcellular location of the 
proteins with differential expression 
(=1.5-fold alterations) in Vero E6 
cells infected with SARS-CoV. A, up- 
regulated proteins. B, down-regulated 
proteins. 


aia cips dase Nuclear, 31, 24.6% 


Unknown location, 


37, 29.4% 
Cytoplasmic, 11, 
Lysosomal, 1, 0.89 8.7% 
Nuclear and 
“cytoplasmic, 4, 
3.2% 


Vesicle, 
Cytoplasmic and~ | | 
Golgi, 1, 0.8% \ 
Vesicle, 1, oxegl | | | Cytoskeletal, 11, 


; 8.7% 
Endoplasmic | 


reticulum, 2, 1.6% | Mitochondrial, 3, 


2.4% 


| 
| 
| 


Secreted or 
extracellular matrix, 
6, 4.8% 


Mitochondrial, 
| peroxisomal and 


. cytoplasmic, 1, 
Membrane protein, | 0.8% 
6, 4.8% 


Down-regulated proteins in SARS-CoV infected cells with at least 


Unknown location, 24, 
40.0% — Cytoskeletal, 5, 8.3% 


hnRNPs are described as a major group of nuclear RNA- 
binding proteins that function in transcription, RNA process- 
ing, MRNA translation, and turnover (21, 22). In the present 


1.5 fold quantitative alterations 


r-Nuclear, 3, 5.0% 
| F Cytoplasmic, 3, 5.0% 
| | Nuclear and 

} / cytoplasmic, 3, 5.0% 


Mitochondrial, 5, 


Vesicular, lyst , 8.3% 


and — 1, \ \ Mitochondrial and 
. ‘o 


\ cytoplasmic, 1, 1.7% 


| 
\_ Membrane protein, 3, 
5.0% 


Vesicle and 
membrane, 1, 1.7% 


Endoplasmic 
reticulum, 4, 6.7% 


Secreted or 
“extracellular matrix, 
7, 11.7% 


work, a total of eight kinds of hnRNP factors were identified to 
be up-regulated significantly in SARA-CoV-infected cells (Ta- 
ble Ill and Supplemental Figs. 4-11), suggesting that the virus 
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Fic. 4. Functional classification of 
the proteins with differential expres- 
sion (=1.5-fold alteration) in Vero E6 
cells infected with SARS-CoV. A, up- 
regulated proteins. B, down-regulated 
proteins. 


Up-regulated proteins in SARS-CoV infected cells with at least 1.5 fold 
quantitative alterations 


Metal binding 
protein, 2, 1.6% 


Translation factor, 5, 
4.0% a 
Ribosomal protein, 
11,8.7% | 
HSP, 2, 1.6% 


hnRNP, 11, 8.7%— 


Cytoskeletal, 3, 2.4% 


Chromatin 
assembly /condensat i— 
on, 4, 3.2% 


Other functions, | 1, } 
8.7% 


Novel protein 
(Function unknown),— 


y- Enzyme, 18, 14.3% 
Signal protein, 16, 
S 12.7% 
Actin network, 9, 
7.1% 


\ Immune response, 5, 
4.0% 


Transcrirtion or 
L replication factor, 9, 
7.1% 


Intracellular 


6, 4.8% | trafficking, 4, 3.2% 


Channel 
protein/transportor,— 
5, 4.0% | 


Protein degradat ion, | 
5, 4.0% 


Down-regulated proteins in SARS-CoV infected cells with at least 
1.5 fold quantitative alterations 


Other functions, 9, _ 


15.0% \ Enzyme, 9, 15.0% 


Novel protein 


(Function unknown), — Signal protein, 10, 


9, 15.0% 16.7% 


Channel 
protein/transportor,—™ 

1, 1.7% | Actin network, 6, 
| Immune response, 8, 10.0% 


Protein degradation, | 
; 13.3% 


2, 3.3% 
Intracellular 


trafficking. 


Transcrirtion or 
replication factor, 3,— 
5.0% 


requires the function of hnRNPs. Among these identified were described here for the first time to be involved in coro- 
hnRNPs, hnRNP A1, hnRNP K, and poly(rC)-binding protein navirus infection (Table Ill). 

have been reported previously to participate in positive-strand HnRNP A‘ has been extensively studied for its role in viral 
virus genome replication. Moreover, the rest of the hnRNPs_ RNA replication. Some reports showed that hnRNP Ai could 
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bind the RNA of a mouse coronavirus (mouse hepatitis virus 
(MHV)) at the 3’-end of both plus and minus strands and 
modulate MHV RNA synthesis (21-23). However, a recent 
work argued that hnRNP A1 might not be necessary for MHV 
viral genome replication or transcription in vivo because the 
absence of hnRNP A‘1 in infected cells had no effect on the 
production of infectious MHV (24). In addition, another exper- 
iment revealed that a mouse erythroleukemia cell line, CB3, 
did not express hnRNP A1 but still supported MHV replica- 
tion, whereas hnRNP A2/B1, hnRNP A/B, and hnRNP A3 
could replace hnRNP A‘1 in cellular functions and viral infec- 
tion (25). In the present study, hnRNP A1 has been identified 
and quantitated by ICAT analysis based on doubly charged 
and both ICAT light reagent- and heavy reagent-labeled pep- 
tide K. YHTVNGHNCEVR.K, a typical trypsin-digested peptide 
fragmentized at the lysine or arginine carboxyl end, which is 
marked with the period (Supplemental Fig. 4) and shows 
up-regulated expression (E6:E6-V, 1:1.90) in SARS-CoV-in- 
fected cells (Table Ill). The present work provides a new 
explanation for such controversial results on hnRNP A1, i.e. 
several different hnRNP factors in the infected cells may form 
a functional hnRNP complex participating in viral RNA metab- 
olism in which one hnRNP factor can be substituted by an- 
other without disruption of the function of the hnRNP 
complex. 

The viral genomic RNA of all positive-strand RNA viruses 
need to be translated by recruited host factors (26, 27). A 
recent report showed that the decrease of 60 S ribosome 
protein levels reduced a positive-strand virus (Brome mosaic 
virus)-directed expression in yeast cells (27). The quantita- 
tively proteomic approaches used here revealed that about 
nine kinds of ribosomal proteins, including components of 
both 40 and 60 S ribosomal subunits, were up-regulated 
significantly (Table Ill and Supplemental Figs. 12-20), sug- 
gesting that overall up-regulation of the ribosomal protein 
expression is required for positive-strand RNA virus propaga- 
tion in host cells. 

Translation factors have been well documented to partici- 
pate in the processes of the virus RNA and protein synthesis 
(26, 28). For example, the elongation factors EF-Tu and EF-Ts 
were found to bind tightly to the viral RNA-dependent RNA 
polymerase (29). We also identified the up-regulation of EF-Tu 
expression in SARS-CoV-infected cells (Table III and Supple- 
mental Fig. 21). Some studies showed that eukaryotic trans- 
lation initiation factor elF-4 protein complex is involved in viral 
protein synthesis (26). On the other hand, the present studies 
revealed that the expression of elF-1, elF-5, and elF-5A2 was 
increased in the infected cells, suggesting that these factors 
play a role in the process of SARS-CoV viral translation (Table 
Ill and Supplemental Figs. 22-24). 

In summary, we present the first quantitative proteomic 
work on the cellular responses to SARS-CoV infection, estab- 
lishing so far the most comprehensive differential proteomic 
index for SARS-CoV-infected cells. The identified differential 


profile derived from the infected cellular proteins gives the 
implications for the infectivity and pathogenesis of SARS-CoV 
and provides a valuable resource for diagnosis, drug devel- 
opment, and clinical treatment for SARS. 
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